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GLOSSARY  OF  NOTATION 


s  cr 


reference  configuration 
vectors  in  1R^  based  at  X  S  S 


B  TR  >  x  -  cJ)(X) 


deformation 


u:  B  +1R 


infinitesimal  deformation 


e  =  [7u+(Vu)  ]/2 


strain 


all  deformations  $ 


F  -  D<j) 


deformation  gradient 


T 

C  =  F  F 


transpose  of  F 
Cauchy  Green  tensor 


Stored  energy  function 


S  =  2  3L 
b  3C 

A  =  — 

3F 

r  JS 
'  3C 

c  =  2CL 


first  Piola-Kirchhoff  stress 


second  Piola-Kirchhoff  stress 


elasticity  tensor 
(second)  elasticity  tensor 


classical  elasticity  tensor 


I  or  L 


identity  map  on  1R  or  8 


%  -  (B,t) 


a  (dead)  load 


all  loads  with  total  force  zero 


LfT^BJR  ) 

sym(T^8,T^B) 

S0(3) 


all  linear  maps  of  T^3  to  IR 
symmetric  linear  maps  of  to  T^5  . 

(Q  G  tOR"5^)  I  QTQ  =  r,  det  Q  =  1} 

im3  m5) 


skew  =  so(3) 


symmetric  elements  of  M„ 

o 

skew  symmetric  elements  of 


infinitesimal  rotation  about  the  axis  v 


equilibrated  loads 


r 


m 


k  :  L  -  M_ 

A  =  kW 

j  =  (k j (ker  k)1)*1 
Skew  *  j(skew') 

Sym  =  j  (syi n) 

$ :  C  -  L 
C 

sym 

Si 

F :  L  -*  Skew 
e 

F:  ]R  *  L  Skew 
e 


astatic  load  map 

astatic  load  for  a  load  2, 

non  singular  part  of  k 

skew  viewed  in  load  space 

sym  viewed  in  load  space 

$(4>)  =  (-DIV  P,  P-N) 

{u:  8  ■+3R'5|u(0)  =  0,  7u(0)  €  sym} 

image  of  C  near  I„  under  $ 

6  sym  8 

W  is  the  graph  of  F 

F(X,Jl)  =  F(XZ)/XJ 

Q's  in  S0(3)  that  equilibrate  A, 
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be  light-face  italic. 


of  notation,  vectors  and  tensors 
are  boldface  italic  and  4-tensors 
other  mathematical  symbols  will 
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§1.  \lntroduction 

x? - 

The  purpose  of  this  paper  is  to  study  the  traction  problem  in 
three  dimensional  nonlinear  elasticity  using  geometric  techniques  and 
singularity  theory.  The  first  two  napers  in  the  series  deal  with  dead 
loads  and  configurations  that  are  nearly  stress  free.  As  was  shown  by 
Signorini  [1950]  and  Stoppelli  [1958],  this  problem  has  nontrivial 
solutions.  However,  their  analysis  is  incomplete  for  three  reasons. 
First,  their  load  was  varied  only  by  a  scalar  factor;  in  a  full 
neighborhood  in  load  space  of  a  load  which  has  an  axis  of  equilibrium 
there  are  additional  solutions  missed  by  their  analysis.  Second, 
their  analysis  is  only  local  in  the  rotation  group,  so  additional 
nearly  stress  free  solutions  are  missed  by  restricting  to  rotations 
near  the  identity.  Third,  some  classes  of  loads  with  a  degenerate  axis 
of  equilibrium  were  not  considered.  This  paper  completes  their  analysis 
by  treating  these  questions  as  well  as  stability.  The  complexity  of 
the  answer  is  indicated  by  the  fact  that  near  certain  types  of  loads, 
we  find  up  to  40  distinct  solutions  that  are  nearly  stress  free.  Our 
constitutive  hypotheses  on  the  stress  tensor  are  ’generic’;  for  a 
degenerate  stress  tensor  there  can  be  even  more  solutions. 

The  literature  on  this  problem  is  very  extensive,  going  back 
to  Signorini  in  the  1930's.  Our  primary  sources  have  been  Stoppelli 
[1958],  Grioli  [1962],  Truesdell  and  Noll  [1965],  Van  Buren  [19t>8], 

Wang  and  Truesdell  [1975],  and  Capri:  and  Podio  Guidugli  [1974], 

However,  none  of  the  references  beyond  Stoppelli  [1958]  proves  any 
of  the  theorems  dealing  with  nontrivial  cases;  i.e.  loads  with  axes 
of  equilibrium.  However,  Grioli  [1962]  is  a  convenient  reference 


for  the  statements. 


The  outline  of  this  first  part  of  the  paper  is  as  follows.  Our 
notation  for  nonlinear  elasticity  and  the  problem  near  a  natural  state 
is  formulated  in  Section  2.  In  Section  3  the  basic  properties  of  the 
astatic  load  are  reviewed  and  developed.  The  problem  is  reformulated 
with  special  reference  to  the  global  aspects  of  the  rotation  group  in 
Section  4  and  introduces  the  bifurcation  equation  which  plays  a  crucial 
role  throughout  the  paper.  Section  5  treats  loads  with  no  axis  of 
equilibrium;  there  are  three  new  features  in  this  section.  First  the 
proof  of  Stoppelli’s  results  is  considerably  simplified.  Second,  the 
results  are  global  relative  to  the  rotation  group.  Finally,  the  stabi¬ 
lity  of  the  solutions  is  determined.  The  number  of  solutions  will  be 
classified  by  load  types;  this  classification  scheme  is  explained  in 
Section  6.  (Some  work  related  to  the  "type  classification”  was  given  by 
Ogden  [1977]).  In  Section  7  a  second  order  bifurcation  equation  is  showai 
to  be  a  gradient.  This  consequence  of  Betti  reciprocity  is  basic  to  our 
analysis.  Section  8  gives  a  complete  bifurcation  analysis  of  type  1 
loads  (the  case  considered  by  Stoppelli),  including  a  stability  analysis. 
New  local  and  global  solutions  are  found.  The  final  section  makes 
explicit  the  comparison  with  Stoppelli’s  theorem. 

The  second  paper  of  this  series  will  analyze  the  remaining  types 
2,  3  and  4,  using  a  reformulation  of  our  gradient  results,  discuss 
linearization  stability,  parallel  loads  and  will  give  additional 

connections  with  the  literature. 
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Statement  of  the  Problem 


b  * 

Let  8  C  TR  be  an  open  bounded  set  with  smooth  boundary  and 

assume  for  convenience  that  ()  E  8.  Let  1  <  p  <  »,  s  >  (3/p)  +  1 

^  5  so 

and  let  C  be  the  space  of  maps  cp:  8  IR  that  are  of  class  W  ** 

1  s  D 

(so  they  are  C  )  such  that  6(0)  =  0,  o  is  a  W  ’^-diffeomor sm 

onto  its  image,  and  J($)  >  0,  where  J  (f )  is  the  Jacobian  of  :< . 

For  example,  if  y:  8  ^  IRJ  is  WS’^-close  to  the  identity  and 

MO)  =  0,  then  ^  G  C  .  If  Q  is  a  linear  isomorphism  of  IRJ  to 

R°  with  det  Q  >  0,  then  Qoii;  EC  as  well. 

Let  points  in  §  be  denoted  X  E  8  and  points  in  be 

denoted  x.  Sometimes  we  write  x  =  o(X).  Let  T^8  be  the  tangent 

space  to  8  at  X,  regarded  as  vectors  in  IRJ  based  at  X.  We  do 

not  identify  T^S  and  IR^  for  conceptual  clarity.  For  c  E  C  ,  let 

F(X)  E  L(T^8,  JRJ)  be  the  derivative  of  6  at  X;  by  standard  abuse 

of  notation  we  write  F(X)  =  D6  (X)  or  76  (X)  interchangeably  and 


L(T^8,  1RJ)  denotes  the  set 

of  all 

linear  maps  of 

TXS 

to 

We  let  F(X)T£  LOR'3 ,  Tx31 

denote 

the  adjoint  of 

F  (X) 

relative  to 

the  Euclidean  inner  product.  Observe  that  F(X)  E  L+(T^8,  IRJ) ,  the 
linear  transformat  ions  with  positive  determinant,  since 
det  F(X)  =  J(3)(X)  >  0.  We  let  C  =  FTF  (that  is, 

CfX)  =  F(X;  F(X)  €  L(T^B,  Tx8)  )  denote  the  Cauchy-Green  tensor. 


We  believe  that  our  results  also  hold  when  8  has  piecewise  smooth 
boundary.  This  program  depends  on  elliptic  regularity  for  such  regions. 
Except  in  special  cases,  this  theory  is  non-existent  and  seems  to 
depend  on  a  modification  of  the  usual  Sobolev  spaces  near  corners. 
However,  for  simple  shapes  like  cubes,  the  necessary  regularity  can 
be  checked  by  hand  in  situations  where  the  linearised  elastostatic 
equations  can  be  solved  explicitly. 


4 


Observe  that  C(X)  £  Svm  (TVB,  TVB)  ,  the  positive  definite 
^  '  pos  X  X  r 

symmetric  linear  transformations  on  T^S  . 

Assume  we  are  given  a  smooth  stored  energy  function  W  defined  on 

Dairs  ( X ,  C )  where  C  G  svm  fTvB,  TVS).  For  ^  -  C  .  the  stored  energy 

'  pos  X  X  - *- 

of  i)  is  |  .W(X,C(X))dVfX) ,  where  C  is  the  Cauchv-Oreen  tensor  of  D 
and  dV  is  the  volume  element  in  5  .  The  fact  that  V.'  depends  on  !> 
only  through- the  ooint  values  of  C  is  usually  called  material  frame 
indifference.  (See  Truesdell  and  Noll  [1965]  and  Marsden  and  Hughes 

[1978]  for  discussions.)  Since  C  is  a  function  of  F,  we  shall, 

T 

by  abuse  of  notation  also  write  W(X,F)  for  W(X,F  F) . 

The  first  Piola-Kirchhof f  stress  tensor  P(X,F)  is  defined 

3 

by  p(X,F)  =  W(X,F) ,  the  partial  derivative  of  W  with  respect 

o  r 

to  F.  Thus,  P(X,F)  £  L(T^SJR^)  .  The  second  Piola-Kirchhoff  stress 

tensor  S(X,C)  is  defined  by  S(X,C)  =  2“^V(X,C),  so 
-  <j  c 

_  * 

S(X,C)  €  sym(T^8,  T^B)  .  From  the  chain  rule,  one  has  the  relationship 
P(X,F)*G  =  4s(X,C)*FTG  +  gtf; 
for  all  G  S  L(T„8  JR'5)  .  7 

For  finite  dimensional  vector  spaces  V,W,  the  bilinear  map 
L(IV,V)  x  L(V,W)  -*•  IR;  (A,B)  *+  ~trace(A  B)  defines  an  isomorphism 

k 

L(V,(V)  =  L(W,V).  If  V  *  W  is  an  inner  product  space,  this 

k 

isomorphism  identifies  sym(V,V)  with  syn(Y,Vl.  Using  this  identifi- 
cation,  we  get  PfX,F)  €  L(1R5,  T^S)  ,  S(X,C)  €  sym(T^S,T^81  and 
PfX.F)  =  S(X,C)oFT,  or  P  =  SFT  for  short, 

"TTie  mass  density  does  not  appear  in  our  formulas  as  we  are  building 
it  into  the  definitions  and  work,  for  example,  with  body  force 
per  unit  volume  rather  than  per  unit  mass. 


D 


Let  A(X,F)  =  denote  the  elasticity  tensor.  Identifying 

0  t*  —  —  ■  -  ■  - 

LCHT^B  JR'21)  ,  L(T^B,  IRJ)*)  in  the  usual  way  for  second  derivatives 
with  B(l(T^6,  IR*5)),  the  space  of  bilinear  maps  of 
L(T^5,  IR0)  *  L{T^8,  IR^)  to  JR  ,  A(X,F)  determines  a  symmetric 
bilinear  map  on  L(T^8,  1RJ)  . 

The  second  elasticity  tensor  C(X,C)  is  similarly  defined  to 


$S 

3C 


be  *  2  evaluated  at  (X,C)  ,  so  may  be  regarded  as  a 


symmetric  bilinear  map  on  >ym(T^5,  T^B)  .  The  chain  rule  gives 


2A(X,F) *(G,H)  =  C(X,C1-(FTH  *  HTF,  FTG  +  GTF)  *  S(X,C)*(HTG  +  GTH) 


The  following  two  assumptions  will  be  made  in  the  first  two 
papers  of  this  series. 

(HI)  The  undeformed  state  is  stress  free;  i.e.  P(X,I)  =  0, 
or  equivalently,  S(X,I)  =  0,  where  I  is  the  identity. 

(H2)  Strong  ellipticity  holds:  there  is  an  e  >  0  such  that 

A(.C,  I)  •  ( v  ®  Z ,  V  ®  Z)  ^  £ll£ii~  II v II “ 

for  all  £  <E  T^B  and  v  €3R*\  where  v  0  $  €  L(T  8,  IR'5)  is  defined  by 

‘  '(v  0 

The  classical  elasticity  tensor  is  defined  by  c(X)  =  -C(X,I)  , 
so  c(X)  is  a  symmetric  bilinear  mapping  on  sym(T^B,T^8)  to  ]R;  at 
$  =  I,  we  identify  T^8  and  1R*5  since  x  and  X  coincide.  By 


A(X, I) • (G,H)  =  jC(X) • (G  +  GT,  H  +  HT) 


(HI)  , 


Regarding  A(X,I)  L(L(T^8,T^8) ,  L(T^8,  T^S))  an^ 
C(X)  S  L(sym(T  8,  TJ),  sym(T  8,  T  6))  ,  this  reads 


2A(X,I)*G  «  c(X)*(G  «•  G  ) 


or,  if  G  is  symmetric, 


A (X , I) *G  =  c(X)(G)  • 


By  (H2),  solvability  of  the  linearized  equations  of  elastostat ics  can 
be  determined  by  the  Fredholm  alternative  (see,  e.g.,  Marsden  and 
Hughes  [ 197S ] ) . 

We  shall  let  B:  8  *►  ^  denote  a  given  body  force  (per  unit 
volume)  and  r:  36  - a  given  surface  traction  (per  unit  area). 
These  are  dead  loads;  in  other  words,  the  equilibrium  equations 
for  $  that  we  are  studying  are: 


DIV  P(X,  F[X))  +  B(X)  =0  for  XE  8 


P(X,F (X) ) • N (X)  -  T(X) 


for  X  E  38 


where  N(X)  is  the  outward  unit  normal  to  38  at  X  E  36  and  DIV  F 
is  the  divergence*  of  P(X,F(X))  with  respect  to  X. 


class  W 


Let  L  denote  the  space  of  all  pairs  (B,t)  =  i  of  loads  (of 
W*  on  5  and  VIs  *  on  38  )  such  that 


B(X)dV(X)  +  !  t ( X ) d A ( X )  =  0 
8  '  38 


*Rccal 1  from  above  that  P(X,F(X))  €  L(R3,  TXS)  .  For  any  v  €  R3 
then  P(X,F(X))v  defines  a  vector  field  on  8  ;  its  divergence 
defines  the  vector  field  DIV  P  hy  (DIV  Pl*v  =  n TV  (PV) 


i . e .  the  total  force  on  8  vanishes,  where  dV  and  dA  are  the 


respective  volume  and  area  elements  on  8  and  3 8.  Using  the 
divergence  theorem  ,  observe  that  if  the  pair  (B,t)  is  such  that  (E) 
holds  for  some  p  G  C,  then  (B,x)  G  L  . 

Throughout  the  paper,  the  group  S0(3)  =  {Q  6  L(IR^JR^)[  Q^Q  =  I 
and  det  Q  =  + 1 }  of  proper  orthogonal  transformations  will  play  a 
key  role. 

By  (HI)  ,  =  Ig  (the  identity  map  on  8)  solves  (E)  with 

B  =  t  =  0.  By  material  frame  indifference,  <p  *  Q j  8  is  also  a 
solution  for  any  Q  G  SO(3) .  The  map  Q  w-  Q 1 8  embeds  SO(3)  into 
C  and  we  shall  identify  its  image  with  SO(3) .  Thus,  the  "trivial" 
solutions  of  (E)  are  elements  of  S0(5). 

Our  basic  problem  is  as  follows: 

(PI)  Describe  the  set  of  all  solutions  of  (E)  near  the  trivial 
solutions  SO(3)  for  various  loads  i  G  i  near  zero. 

Here  "describe”  includes  the  following  objectives: 

(a)  counting  the  solutions 

(b)  determining  the  stability  of  the  solutions 

(c)  showing  that  the  results  are  insensitive  to  small 
perturbations  of  the  stored  energy  function  and  the  load. 
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S3 .  The  Astatic  Load  and  .Axes  of  Equilibrium 

This  section  is  devoted  to  the  geometry  of  the  load  space  L . 
Many  of  the  results  of  this  section  are  available  in  the  literature, 
but  we  gather  them  here  for  convenience. 

Before  beginning,  we  shall  recall  a  few  notations  and  facts 
about  the  rotation  group  SO ( 3 j .  Let 

=  KIR'5,  IR5)  =  linear  transformations  of  IR5  to 
sym  =  {A  S  M,|AT  =  A} 
skew  =  {A  S  M„|A^  »  -A} 


We  identify  skew  with  so(3),  the  Lie  algebra  of  S0(3);  skew  and 

3  3  /V 

F  are  isomorphic  by  the  mapping  v  E1R  v  E  skew  ,  where 

v(w)  =  w  x  v.  If  v  =  (p,q,r)  relative  to  the  standard  basis,  then 


f  0  r 


v 


-r 


Lq 


0 

-P 


-q 

p 

o 


The  Lie  bracket  is  [v , w ]  =  v  ®  w  -  w  ®  v  =  - (v  *  w)  where 

v  K  w  G  M,  is  given  by  (v  X*  w)  (u)  =  v<w,u)  .  The  inner  product 
1  N  t 

is  <v,w>  =  —trace  (v  w)  ,  the  Killing  form  on  so(5).  Finally, 
exp(v)  is  the  rotation  about  the  vector  v  in  the  positive  sense, 
through  the  angle  II v It  - 

Now  we  turn  to  a  study  of  L  .  For  t  €C  and  l  E  L  ,  we  say 
that  l  is  equilibrated  relative  to  p  if  the  total  torque  in  the 
configuration  <p  vanishes: 


[  <*>(X)  x  BlX)dV(X)  +  [  (J)fXJ  x  T(X)dA(X)  =  0 
JS  ’  38 


where  i  =  From  symmetry  of  the  stress  one  sees  that  if 

l  =  (B,t)  £  L  satisfies  (E)  for  some  {>  €  C,  then  Z  is  equilibrated 
relative  to  >J>.  (An  easy  proof  uses  the  Piola  transform;  cf.  Marsden 
and  Hughes  [ 1978] ) . 

Let  L denote  the  loads  that  are  equilibrated  relative  to  the 
identity  configuration  Ig  . 

Define  the  astatic  load  map  k:  L  *  C  by 


k(Z,<p) 


B(X)  X>^>(X)dV(X)  ♦  [  T (X)  ».J>(X)dA(X) 

JS  j3B 


and  write  k(£)  »  k(Jl,  Ig) . 

We  have  actions  of  SO (5)  on  L  and  C  given  by: 


Action  of  SO (5)  on  L:  Q£(X)  =  (QB(X) ,  QtfX)) 
Action  of  SO (5)  on  C :  =  Q0(P 


Note  that  QZ  means  ,fthe  load  arrows  are  rotated,  keeping 

the  body  fixed.”  We  shall  write  0^  and  C ^  for  the  SO(5)  orbits 

of  Z  and  <fc.  Thus,  0  denotes  the  trivial  solutions  corresponding 

1  8 


to  Z  =  0. 


The  following  is  3  list  of  basic  observations  about  the  astatic 


load,  each  of  which  is  readily  verified: 
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t Ai )  l  is  equilibrated  relative  to  $  if  and  only  if 

k(x,,4>)  €  sym.  In  particular,  l  E  if  and  only  if 
k(Z)  E  sym. 

(A2)  (equivariance)  .  For  l  E  L,  D  €  C,and  ,Qn  E  S0(5) 

k(Qii.Q2-»  =  Q^U.^OqT 

In  particular,  k(Q£)  =  Qk (£) 

(A3)  (infinitesimal  equivariance).  For  £E  L,  d)  E  C  ,  and 
W  E  skew, 

k(W^,f)  =  k (il , =  -kU,$)  W\ 

In  particular,  k (W£)  =  Wk(2). 

Later  on,  we  shall  be  concerned  with  how  the  orbit  of  i  E  L 
meets  i  ,  The  most  basic  result  in  this  direction  is  the  following. 

5.1.  DaSilva's  Theorem.  Let  l  E  i.  Then  ().  n  [  f  0, 

-  -  - _  l  e 

Proof.  By  the  polar  decomposition,  we  can  write  k(/l) 

for  some  QE  S0(5)  and  AE  sym.  By  (A2)  ,  k(Q£)  =  Qkl<A 

so  by  C A 1 ) ,  Q£  E  L  .  ■ 
e 

Similarly,  any  load  can  be  equilibrated  relative  to  any  chosen 
configuration  by  a  suitable  rotation. 

The  concept  of  an  axis  of  equilibrium  deals  with  the  case  in 
which  C ■  meets  i ^  in  a  degenerate  way. 


T 

*  Q  A 


=  A  E  sym, 
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5.2.  Definition.  Let  l  G  L  and  v  GIR'5,  llvll  =  1.  We  sav  that 
-  e 

v  is  an  axis  of  equilibrium  for  Z  when  exp(6v)£  G  L for  all 
real  9,  i.e.  when  rotations  of  £  about  the  axis  v  do  not 
destroy  equilibration  relative  to  the  identity. 

A  number  of  useful  ways  of  reformulating  the  condition  that  v 
be  an  axis  of  equilibrium  are  as  follows. 

3.5.  Proposition.  Let  Z  G  and  A  =  k(£)  6  sym.  The  following 
conditions  are  equivalent: 

1 .  Z  has  an  axis  of  equilibrium  v 

2.  there  is  a  v  G  1R*5  ,  llvll  =  1  such  that  vZ  G  L 

3.  W  h*  AW  *  WA  fails  to  be  an  isomorphism  of  skew  to  itself 

4.  trace  A  is  an  eigenvalue  of  A. 

Proof.  1  =*  2.  This  follows  by  differentiating  expi5v).\  in  9 
at  9=0.  2  =*  1.  By  (A2) , 

k(exp(9v)?J  =  'I  +  (9v)  +  4{(9v)“  +  ...]k[Z)  , 

Since  k(v’)  =  vk(£)  is  symmetric,  this  is  symmetric,  term 
by  term. 

2  =**  3.  Since  k(v»l)  =  vA  is  symmetric,  vA  +  Av  =  0,  so 
W  K  AW  +  WA  is  not  an  isomorphism. 

5  ^  2.  There  exists  a  v  £  ^  ,  llvll  =  1,  such  that  vA  +  Av  =  d, 

so  k(vo)  =  vA  is  symmetric. 

3  4 .  Define  L  €  by  L  =  (trace  AjI  -  A.  Then  one  has 


the  relationship 
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(Lvr  *  Av  *  vA 

(In  fact,  if  [u,v,w]  denotes  the  triple  product,  the  relationship 
[Bu,Bv,Bw]  =  det  B[u,v,w]  gives  [Au,v,w]  +  [u,Av,  v;]  +  [u,v,Aw] 

=  (trace  A)  [u,v,w]  .  This  yields  (Lv)  =  A\A  +  A  v  ,  which 

gives  the  claimed  results  for  symmetric  A).  Therefore, 

Av  +  vA  =  0  if  and  only  if  Lv  =  0  ,  i.e.,  v  is  an  eigenvector 
of  A  with  eigenvalue  trace  A  .  • 

3.4.  Corollary.  Let  £  €  Lq  and  A  =  k (il)  E  svm.  Let  the  eigen- 
values  of  A  be  denoted  a,  b,  c.  Then  £  has  no  axis  of  equili- 
brium  if  and  only  if  (a  +  b) (a  +  c] (b  +  c)  ^0. 

Proof.  This  condition  is  equivalent  to  saying  that  trace  A 
is  not  an  eigenvalue  of  A.  ® 

5.5.  Definition .  We  shall  say  that  ll  G  [g  is  a  type  0  load  if  l 
has  no  axis  of  equilibrium  and  if  the  eigenvalues  of  A  =  k(£) 
are  distinct. 

The  following  shows  how  the  orbits  of  type  0  loads  meet  L  . 

3.6.  Proposition .  Let  l  E  L&  be  a  type  0  load.  Then  Q  ^ 
consists  of  four  type  0  loads. 

Proof.  We  first  prove  that  the  SO (5) -orbit  of  A  in  under 
the  action  Q~~QA  meets  sym  in  four  points.  Relative  to  its  basis 
of  eigenvectors,  we  can  write  A  *  diag(a,b,c).  Then  0^  ~  ^yn 
contains  the  four  points 


diag(a,b,c) 
diag(-a,-b,c) 
diag(-a,b , -c) 
diag(a , -b , -c) 


(Q  =  I) 

(Q  =  diagC-1,-1,1)) 
(Q  =  diag(-l , 1- 1) ) 
IQ  =  diag(l 1 , - 1) ) 


These  are  distinct  points  since  (a  +  b)  (a  +  c)  (b  +  c)  0.  Now 
suppose  a,  b  and  c  are  distinct.  Suppose  QA  =  S  G  sym.  Then 

S“  =  A-.  Let  be  an  eigenvalue  of  S  with  eigenvector  u. . 

2  9  2  °  n 
Then  S  =  y7u^  ~  so  an  eigenvalue  of  A“ .  Thus, 

2 

as  the  eigenvectors  of  A“  with  a  given  eigenvalue  are  unique, 
u^  is  an  eigenvector  of  A  and  is  the  corresponding  eigenvalue. 

Since  det  Q  «  +1,  det  S  -  det  A,  so  we  must  have  one  of  the  four 
cases  above. 

By  equivariance ,  kCc)^)  ^  sYm  ~  ^k(£)  ^  S^TT1  a  set  cons^st^-ng  °f 
four  points.  Now  0 ^  "  L0  =  k  *(0^^  ^  svm) ,  so  it  suffices  to 
show  that  k  is  one-to-one  on  0^.  This  is  a  consequence  of  the 
following  and  (A2) . 


3.7.  Lemma .  Suppose  A  £  sym  and  dim  ker  A  ^  1.  Then  A  has  no 
isotropy :  i . e,  QA  -  A  implies  Q  =  I. 

Proof.  Every  Q  i  I  acts  on  IR*5  by  rotation  through  an  angle 
9  about  a  unique  axis,  that  is,  about  a  line  through  the 
origin  in  IR^  .  Now  QA  =  A  means  that  Q  is  the  identity  on  the 
range  of  A.  Therefore  if  0  i  I  and  QA  =  A,  the  range  of  A 
must  be  0  or  1  dimensional,  i.e.  dim  ker  A  7*  2.  ■ 
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Finally  in  this  section,  we  study  the  range  and  kernel  of 
k:  L  -*•  M- 

o 

3.8.  Proposition.  1.  ker  k  consists  of  those  loads  in  L for 
which  every  axis  is  an  axis  of  equilibrium. 

2.  k:  L  +  M-  is  surjective. 

J  - i - 


Proof.  1.  Let  i  G  ker  k.  For  W  €  skew,  k(W&)  =  Wk(£)  =  0 

so  W£  G  Lq;  by  3.3  every  axis  is  an  axis  of  equilibrium.  Conversely, 

if  €  l  for  all  W  G  skew,  then  k(W£0  =  Wk(20  is  svmmetric 
e 

for  all  W;  i.e.  k(Ji)W  +  Wk(Z)  =  0  for  all  W.  From 

(Lv)^  =  Av  +  vA,  where  A  =  k(£) ,  and  L  =  (trace  A) I  -  A,  we 

see  that  L  =  0.  This  implies  trace  A  =  0  and  hence  A  =  0. 

To  prove  2  introduce  the  SG(3)  invariant  inner  product  on  L: 


<  1,1) 


r  f 

(  B(X)  ,B(X)  >dV(X)  +  |  (  t(X)  ,t(X)  >dAl.X)  . 

J3B 


T 

Relative  to  this  and  the  inner  product  trace  (A  Bj/C  on  NL,  the 
T 

adjoint  k  :  M  L  of  k  is  given  by 


kT( A)  =  (B,t)  where  B(Xj  =  AX-G,  tfX)  =  AX-G, 

and 

G  =  ft  AX  UV(X)  +  f  AX  dA(X)1/: jdV  *  '  JA  ^  . 

1  '-B  J  “  8  •' 3  8  J 

If  kT(A)  =  (0,0)  then  it  is  clear  that  A  =  0.  It  follows  that 
k  is  surjective.  8 


I 
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5.9.  Corollary .  1.  ker  k  is  the  largest  subspace  of  that  is 

S0(3)  invariant  . 

2.  k]  (ker  k)^  :  (ker  k)^  ->■  is  an  isomorphism. 

Let  j  =  (k|(ker  k)^)~^  and  write 

Skew  =  j  (skew)  Sym  =  j  (sym) 


These  are  linear  subspaces  of  L  of  dimension  3  and  6  respectively. 
Thus  we  have  the  decomposition: 

S0(3) 

1  _ 

L  =  Skew  0  Sym  ®  ker  k 

L 


corresponding  to  the  decomposition  NL  =  skew  <£  sym; 

1  T  1  T 

u  =  ’  U  )  +  y(U  +  U  ) . 


Note:  Skew  and  L  need  not  be  orthogonal. 

-  e 


J4 .  Equivalent  Reformulations  of  the  Problem 


Define  0:  C  L  by  $(?)  =  (-DIV  P,  Pv\T)  i.e. 


«($)(X)  =  (-DIV  P(X,F(X))  ,  P (X , F (X) )  *N (X) ) 


so  the  equilibrium  equations  {E)  become  $($)  =  l  .  From  material 
frame  indifference  we  have  equivariance  of  4 ;  i.e.  5(Qc)  =  Q-M  $)  . 
Standard  Sobolev  estimates  show  that  $  is  a  smooth  mapping  (see, 
for  example,  Palais  [1968]).  The  derivative  of  b  is  given  by 


D$(cp)*u  s  (-DIV(A* Vu)  ,  (A*7u)  *NJ 


and  at  <p  =  Ig  this  becomes 


»(Ig)*u  =  (-DIV(C*e),(C*e)-N) 


1  T  - 

where  e  =  — [7u  +  (7u) 

I f  D$(Ir)  •  Tt  C  L  were  an  isomorphism  we  could  solve 

5  8 

$($)  =  l  uniquely  for  $  near  Ig  and  l  small.  The  essence 

^  ....  .  .  y since  ^ (SO (3 ) )  =  0, 

of  our  problem  is  that  D<P(Ig)  is  not  an  isomorphism;  kernel’ D*(TT"ccTritain 

Define  C  =  fu  £  Tt  C1  u(0)  =  0  and  7u(0)  £  svm}  .  From 
sym  Ig  ' 

(H2)  and  linear  elasticity,  we  have: 


4.1.  Lemma.  D$fIJ  'C  :  C 

-  8  '  svm  sym 


Le  is  an  isomorphism. 


I 


The  connection  between  the  astatic  load  map  k:  L  -*■  M  and  * 
is  seen  from  the  following  computation  of 

-1  2.  Lemma.  Let  v  G  C  and  let  P  be  the  first  Piola-Kirchnoff 
stress  tensor  of  .  Then 

kf }(':))  =  '  ?dV 

3 

This  follows  by  an  application  of  Gauss'  theorem  to 

f  f 

k($0'n  =  |  C-DIV  P)  ®  X  dV(X)  +  '  fP  \-')  »  X  dA(Xl. 

'  8  j36 

This  should  be  compared  with  the  astatic  load  relative  to  the 
configuration  b  rather  than  I„ ;  one  gets 

k($($)  ,5)  =  !  S  dV 
;S 

which  is  symmetric,  while  k(*>(£))  =  k f$ (40  ,  I_)  need  not  be. 

o 

To  study  solutions  of  $(4>)  =  i  for  4>  near  the  trivial 

solutions  and  ?,  near  a  given  load  2^,  it  suffices  to  take 

■'q  ~  L  .  This  follows  from  Da  Silva's  theorem  and  equi variance  of 

Let  Csvm  be  regarded  as  an  affine  subspace  of  C  centered  at 

In.  Let  ?  be  the  restriction  of  £  to  C  .  From  the  implicit 
8  sym  v 


function  theorem  we  get: 


Later  we  shall  show  how  to  compute  D“F(0)  in  terms  of  0$(lg) 
and  C. 

Now  we  are  ready  to  reformulate  problem  (PI). 

(P2)  For  a  given  l ^  E  near  zero,  study  how  0 ^  meets  the  graph 
of  F  for  various  l  near  ^ . 

Problems  (PI)  and  (P2)  are  related  as  follows.  Let  o  solve  (C) 
with  l  E  JL  and  Q  be  such  that  J  =  Q<j>  E  C  .  Then  $(0)  =  Q£  , 

x  xv  sy7n 

so  the  orbit  of  i  meets  the  graph  of  F  at  £($).  Conversely,  if 
the  orbit  of  l  meets  M  at  $(5)  =  Q£  >  then  (j>  =  Q  solves  (E). 
We  claim  that  near  the  trivial  solutions,  the  numbers  of  solutions 
to  each  problem  also  correspond.  This  follows  from  the  next  lemma. 


4,4.  Lemma.  There  is  a  neighborhood  U  of  Ig  in  C such  that 
0  E  U  and  Q<j>  6  U  implies  Q  =  I g . 

Proof.  Note  that  C _  is  transverse  to  CT  at  Ig  and  IR 


sym 


I, 


has  trivial  isotropy.  Thus,  as  S0(5)  is  compact,  0  is  closed,  so 

3 

there  is  a  neighborhood  of  of  Ig  in  ^vm  such  that 

Qj  8  E  Un  implies  Q  =  I.  The  same  thing  is  true  of  orbits  passing 


through  a  small  neighborhood  of  I  by  openness  of  t ransversality  and 

8 

compactness  of  S0(3) .  ■ 

If  0  meets  M  in  k  points  =■  4>(uk),  i  =  l,...,k  then 

5.  are  distinct  as  $  is  1-1  on  a  neighborhood  of  IL,  in  C 
i  b  sym 

If  this  neighborhood  is  also  contained  in  U  of  4.4,  then  the  points 


Q i  5 1  =  ^  are  also  distinct  by  -1.4.  Hence  the  problems  (PI')  and  ;PJ) 

are  equivalent. 

In  connection  with  the  action  (Q, A)  k  QA  of  S0(3)  on  M,  we 
shall  require  some  more  notation.  Let 

Skew  (A)  =  *  A^l  £  skew 

and 

Sym  (A)  =  +  A^  1  E  sym 

be  the  skew  symmetric  and  symmetric  parts  of  A,  respectively. 

We  shall,  by  abuse  of  notation,  suppress  j  and  identify  Sym 
with  sym  and  Skew  with  skew.  Thus  we  will  write  a  load  l  E  L  as 
i  =  CA,n)  where  A  *  k(£)  E  M_  and  n  E  ker  k;  hence  Z  E  precisely 
when  A  E  sym.  The  action  of  SO  (3)  on  L  is  given  by 

Qe  -  (QA, Qn). 

Using  this  notation  we  can  reformulate  problem  (P2)  as  follows: 

(P3)  For  a  given  l n  =  (AiA,nrt)  E  L ,  near  zero,  and  l  =  (A,n) 

-  -  —  - * —  ()  i)  l)  e  — — — -  ■ 

near  find  Q  E  SO (31  such  that 

Skew(QA  F (Sym(QA) ,Qnl  =  0. 

Define  the  rescaled  map  P :  1R  *  Skew  by 


(3.2a) 


(3.2b) 
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F(A,2)  =  — tF l \ 2 ) . 

Since  F(0)  =  0  and  DF(O)  =0,  F  is  smooth.  Moreover,  if 
F(£)  =  4g(&)  +  ^C(Z)  +  ...  is  the  Taylor  expansion  of  F  about 
zero,  then  F(\,£)  =  +  ^C(c)  +  ...  . 

In  problem  (E)  let  us  measure  the  size  of  l  by  the  parameter 
\  .  Thus,  replace  4> (<i> )  =  l  for  Z  near  zero  by  5(0)  =  \l 
for  A  near  zero.  This  scaling  enables  us  to  conveniently  distinguish 
the  size  of  i  from  its  'orientation'.  In  the  literature  l  has 
always  been  fixed  and  A  taken  small.  Here  we  allow  l  to  vary  as 
well.  Thus  we  arrive  at  the  final  formulation  of  the  problem. 

(P4)  For  a  given  l ^  =  (A^,n^)  <=  L  ,  for  l  near  and  \  small , 

find  Q  £  S0(5)  such  that 

Skew(QA)  -  A F  ( \ ,  Sym(QA) ,Qn)  =  0. 

The  left  hand  side  of  this  equation  will  be  denoted  H(\,A,n;Q) 


or  H(A,Q)  if  A,  n  are  fixed. 


-k 


§5 ,  Type  Q:  No  axis  of  equilibrium 

We  shall  begin  the  analysis  by  giving  an  (almost  trivial)  proof 

ft 

of  one  of  the  basic  theorems  of  Stoppelli  ’ 1938' 

5.1.  Theorem.  Suppose  l  E  L  has  no  axis  of  equilibrium.  Then  for 
a  sufficiently  small,  there  is  a  unique  q  E  ^svm  and  a  unique  Q 
in  a  neighborhood  of  the  identity  in  S0(5)  such  that  =  Q  *  f  solves 
the  traction  problem 


P((p)  *  U  . 

Proof.  Define  H:  1R  x  S0(3)  ■+•  Skew  by 

H(X,Q)  =  Skew (QA)  -  XF (X ,Sym(QA) ,  Qn) ) 

where  l  =  (A,n)  E  -  Sym  Q  ker  k  is  fixed.  Note  that 
D^H(0 , 1)  •  W  =  Skew  (  WAJ  =  ^-(WA  +  AW).  By  Proposition  3.5,  this  is 

an  isomorphism.  Hence,  by  the  implicit  function  theorem,  H(X,Q)  =  0 
can  be  uniquely  solved  for  Q  near  I  E  S0(3)  as  a  function  of  X 
near  0  €  IR,  ■ 

The  geometric  reason  "why"  this  proof  works  and  the  clue  to 
treating  other  cases  is  the  following. 


The  only  other  comp lete  proof  in  English  we  know  of  is  given  in 
Van  Buren  tl968’t  although  sketches  are  available  in  (Irioli  19t>3  .  , 
Tmesdell  and  Noll  "19b5^and  Wang  and  TruesJell  .  Our  proof 

is  rather  different;  the  use  of  the  map  F  avoids  a  series  of 
complicated  estimates  used  hv  Stoppelli  and  Van  Buren. 


D  .  - 


Lemma.  A  load 


l  €  has  no  axis  of  equilibrium  precisely 
when  L  -  $  T ^0 ,  .  In  particular,  if  l  has  no  axis  of  equili¬ 

brium.  then  C,  intersects  L  transversely  at  l 
-  I - e  - 1 - 

Proof.  The  tangent  space  to  0 ^  at  l  £  is  T ^0 *  =  WJ1MV  -  skew; 
and  the  projection  of  this  into  the  complement  Skew  to  is 

W2,  *+  h(WA  *  AW)  where  A  =  k  (i)  .  The  result  then  follows  from 
part  3  of  3.3.  • 


We  have  shown  that  there  is  only  one  solution  to  $($)  =  \l  near 
the  identity  if  A  is  small  and  l  has  no  axis  of  equilibrium.  How 
many  solutions  are  there  near  the  trivial  solutions  S0(3)?  As  we 
shall  see,  this  problem  has  a  non-trivia!  answer  which  depends  on  the 

type  of  i.  We  analyte  the  simplest  case  here.  Recall  from 

definition  3.5  that  a  load  'l  E  L  is  said  to  be  of  type  0  if  l  has 
no  axis  of  equilibrium  and  if  A  *  k[l)  has  distinct  eigenvalues. 

Loads  with  no  axis  of  equilibrium  occur  amongst  other  types  of 
loads  classified  in  the  next  section,  and  Stoppelli’s  theorem  3.1 
applies  to  them.  However,  the  global  structure  of  the  solutions 
(‘’'global”  being  relative  to  S0(3))  is  different  for  the  different  types. 
For  type  0  the  situation  is  as  follows. 

5.3.  Theorem.  Let  £  i  be  of  type  0.  Then  for  >  suffic ient  ly 

sma  1 1  ,  $  f  $ )  =  A  ?  q  has  exactly  four  solutions  r  ^  •  t  *  t  -  and  ;•  ^  i  n  a 


neighborhood  of  the  trivial  solutions  SO (3)  C  C 


(  see  Figure 


Proof.  By  3.6,  0,  ,  meets  L  in  four  points.  By  5.1,  in 

vx,o  e 

a  neighborhood  of  0  in  L,  0-.  ,  meets  V  in  exactly  four  points, 

AJO0 

the  images  of  5^,  iL  and  say.  Thus  problem  (P2)  has  four 

solutions.  By  the  equivalence  of  (PI)  and  (P2) ,  so  does  (PI).  ® 


S0(3) ;  the 
trivial  solutions 


Let  A  =  and  5^  =  { Q j QA  G  svm}.  From  the  proof  of  i.b 

we  see  that  is  a  four  element  subgroup  of  SO ( 3 )  isomorphic  to 

®  Z~.  By  our  earlier  discussions,  the  elements  p.  are  obtained 

-  1 

from  <jL  by  applying  rotations  close  to  elements  of  S^.  In 
particular,  as  \  -+■  0,  the  solutions  {-jk}  converge  to  the  four 
element  set  (regarded  as  a  subset  of  C) . 

For  i  sufficiently  close  to  i  the  problem  M:)  =  \l  will 
also  have  four  solutions.  Indeed  by  the  openness  of  transversal itv , 

0^  will  also  meet  .V  in  four  points.  In  other  words,  the  picture 
for  type  0  in  Figure  2  is  stable  under  small  perturbations  of 
Next  we  study  the  stability  of  the  four  solutions  found  in 
Theorem  5.5.  This  will  be  done  under  the  hypothesis  that  the  classical 
elasticity  tensor  in  stable.  We  introduce  the  following  condition; 

(H3)  Assume  there  is  an  n  >  0  such  that  for  all  e  6  SymlT^S,  T^S), 

l  t  t 

c(e)  =  ^c(X)(e,e)  >  pilell^,  il  •  II  =  pointwise  norm 

(e(e)  is  the  stored  energy  function  for  linearized  elasticity). 

Because  of  difficulties  with  potential  wells  and  dynamical 
stability  in  elasticity  (see  Knops  and  Wilkes  ” 1 9T 2 ]  and  Ball,  Knops 
and  Marsden  [ 19^8]  )  we  shall  adopt  the  following  ’’energy  criterion” 
definition  of  stability. 

5.4.  Definition .  A  solution  u  of  M-f)  =  l  will  be  called  stable 
if  p  is  a  local  minimum  in  C  of  the  potential  function 


2b 


V ,  i  ? )  =  W  i  .3 )  dV  - 

J 

B 

r  r 

where  B(X)  (\)dV  (X)  +  T  IX)  •  3  (X  J  dA(X)  =  trace  k  (  h  ,  t 

>  3  ;  )8 

If  0  15  not  stable,  its  index  is  the  dimension  of  the  largest 

subspace  of  vectors  u  tangent  to  C  at  3  with  the  property 
that  ^  decreases  along  some  curve  tangent  to  u .  (Thus,  index  0 
corresponds  to  stability.) 

5.5.  Theorem .  Assume  (Hl)-fH5)  and  let  l ^  he  as  in  5.5.  For 

sufficiently  small,  amongst  the  four  solutions  t-.,  2 ,  given 

— - — - - — - - -  1  _  b  4  - 

by  5.5,  exactly  one  is  stable;  the  others  have  indices  !,  2,  and  5. 
Suppose  i>  is  a  solution  approaching  Q  £  as_  \  -  0.  Then  t 
is  stable  if  and  only  if  QA  -  trace (OA) I  £  svm  is  positive-definite. 
In  general,  the  index  of  $  is  the  number  of  negative  eigenvalues 
of  QA  -  trace  (QA) I  . 

Proof.  Let  t  C  solve  Me)  =  .  Then  is  a  critical 

point  of  V.  .  Consider  the  orbit  c\  =  ;Qv0  Q  -  SOfoV  of  ;  . 
w'0  *0 
Its  tangent  space  decomposes  T,  C  as  follows: 

•:o 

Tv  C  =  t\  •£  IT.  0.  'x 
*  0  -0  *n  "o  *n 

First  consider  V.  restricted  to  fT#  0 t  ) 

*  *0  Ml 

at  in  the  direction  of  a  -  i?  2,  ■ 

0  "  1)  M 


second  ler: .at 


-i 


■  u , 


cfX) ♦ (e(X) ,e(X) jdV(X) ,  where 

8 

1  T 

e  =  —  (7u  +  (7u)  ).  This  is  larger  than  a  positive  constant  times  the 

7 

L~  norm  of  e,  by  (H3)  .  However,  since  u  is  in  (T  C\)% 

2  2 

He  11%  >  (constant)  II  u  ir  .  by  Korn's  ineaualitv  (see  Fichera  1972*). 

L“  H 

By  continuity,  we  have 


p  =  Q|g  ,  this  becomes 


D“V,  ,  (p  )-(u,u)  ^5  Hull 

Axa  U 


if  u  is  orthogonal  to  0  at  and  \  is  small.  This  implies 

*0  0 


that 


'0 


is  a  minimum  for  V, ,  in  directions  transverse  to  C 


"0 


'0 


(Actually  one  can  see  that  is  a  local  minimum  in  the  topology 

i 

of  C  on  (T  0  )“*  bv  using  Tromha’s  "2976  version  of  the  Morse 

1)  !)  * 

lemma. ) 


Next,  consider  V. .  restricted  to  C  .  Bv  material  frame 

» 2  t) 

0  J0 

indifference,  W  is  constant  on  0  and  so  as  must  be  a 

>T  ( ) 


"0 


critical  point  for  V.  restricted  to  0,  ,  it  is  also  a  critical 


point  for  =  l 


0  "0 

estricted  to  0  (where  It)  =  (*',;->!.  It 

p0 


suffices  therefore  to  determine  the  index  of  l  0  at  .  The 


•'  0 


(V 


result  is  now  a  consequence  of  continuity  and  the  limiting  case 
•  -  0  given  in  the  following  lemma  about  type  0  loads. 


5 .  b .  Lemma .  Let  l  he  type  0  and  let  A  =  k I ? ) .  Then  S  , 

regarded  3s  a  subset  of  C  equals  the  set  of  critical  points  of 

0 T  .  These  4  critical  points  are  nondegenerate  with  indices 
8 

f) ,  1 ,  2  and  3 ;  the  index  of  Q  is  the  number  of  negative  eigenvalue^ 


of  >\  - 1 race fOA) I . 


I 

Proof.  First  note  that  L  =  (T,  SO  ( j  1  )“  since  Mil,,)  nas 
*  e  ip  b 

b 

kernel  T  SCH5)  =  skew,  has  range  L  and  is  selt-adi oint .  Thus 
i  *1  e 

o 

Q£  G  L  if  and  onlv  if  l  i  T  S0f5).  It  follows  that  QZ  €  L  if 
e  l  e 

j  0 

and  only  if  Q  is  a  critical  point  of  i’(?T  .  Recall  that 

1  5 

elements  of  3^  =  {Q  E  S0f5):Q2.  6  i  }  are  symmetric.; 

To  compute  the  index  of  </Q  at  QGS  we  comvu te  the 

8 

second  derivative 


dt 


2*  l exp  i'tlVlQ)  1 1=0  =  UiVQ ) 


Now 


v(W“Q)  =  trace  k(M,l\TQ)  -  trace [k f  ’  ,OVkV~ ] 


.  i  “> 

-  trace  [  AQ  W“  ]  5:  trace  [W“QA] 


because  Q  =  Q  .  Tills  Quadratic  form  on  skew  is  represented  bv  the 
element  OA  -  t race (OA )  I  o+  sym  as  is  seen  fror.  vA  +  Av  =  1  Lv  1  *  with 
A  replaced  by  OA  and  trace  ivwl  -  2v*w  .  ,!siny  the  representations 
for  { QA  ’  given  Proposition  S.n,  namely 

diagfa,b,c  ),  di ag * -a , -b  , c  1  ,  di ag ( -a , b , -c )  and  di ag ( a , -n , -c 1 


one  checks  that  all  four  indices  occur.  ■ 


Remark .  This  lemma  is  a  special  case  of  the  general  problem 


to  study  the  critical  points  ;;ne(r  *oinct  j  on".  1  <  orbits  of 


a 


) 


"Hi 


representation  of  a  Lie  group.  This  situation  will  arise  again  in 
our  analysis  of  the  other  load  types,  cf.  Frankel  [1965]  aac 
Ramanujam  ^19691. 


:b.  Class i ficat ion  of  Orbits  in  M.. 

- 0 

The  purpose  of  this  section  is  to  classify  orbits  in  under 

the  action  (Q,A)  QA  of  S0(3J  on  NL  by  the  way  the  orbits 
meet  svm.  It  is  enough  to  consider  orbits  0^  of  elements  of  sym 
by  the  polar  decomposition.  We  begin  by  recalling  a  result  already 
proved. 

b.l.  Proposition  (Type  0).  Suppose  A  G  sym  has  no  axis  of  equilibrium 
and  has  distinct  eigenvalues.  Then  0^  ~  sym  consists  of  four  points 
at  each  of  which  the  intersection  is  transversal. 

This  was  Proposition  5.6.  (Another  proof  of  this  is  given  below.) 

We  shall  let  the  eigenvalues  of  A£  sym  be  denoted  a,b,c. 

Using  the  terminology  from  §5,  we  say  that  A  has  no  axis  of  equilibrium 
when  (a+b) (b  +  c) (a+c)  f  0;  i.e.  a  +  b  +  c  i  a,  b  or  c,  and  in  this 

case  0  intersects  sym  transversely  at  A. 

6,2.  Definition ,  We  shall  say  A  is  of  type  1  if  A  has  no  axis 

of  equilibrium  and  if  exactly  two  of  a,b.c  are  equal  and  non-zero 
i say  a  =  b  ^  c ,  a  ^  0 ) . 

6.5.  Proposition .  I f  A  is  type  1,  then  0 ^  sym  consists  of 

two  points  (each  with  no  axis  of  equilibrium)  and  an  IRT^  ^  ^  ' oa on 

point  of  which  has  one  axis  of  equilibrium) . 

Before  proving  this,  we  give  a  number  of  lemmas  of  general 

__  5 

utility.  If  l  t  1RP"  is  a  line  through  the  origin  in  IK  , 
be  the  rotation  through  angle  t  about 


let  0 
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r. .4.  Lemma.  *  -  Q,  is  an  embedding  of  MP~  onto  S0i3)  ~  sym\  I  . 

Proof.  It  is  clear  that  3  h-  Q,  is  a  one-to-one  nap  of  RP~ 

■»  -IT 

into  S0(5).  Since  Q7  =  I ,  Q,  =  Q,  =  Q.*,  so  it  maps  into 

^  <>, 

SO  (.3)  '■  sym. 

Every  Q  -=  S0i3')\  I  is  a  rotation  through  some  angle  -?  about 
some  axis  <3.  If  such  Q  also  is  symmetric  then  it  has  three 
independent  real  eigenvectors.  Hence  r  =  t  .  ■ 

b.5.  Corollary.  The  orbit  of  the  identity,  (?T,  meets  sym  m 

i 

one  point  (I)  and  FP  =  (S0(3)  sym'A  I. 


6.6,  Lemma .  Let  A  £  sym  with  dim  ker  A  i  1  and  suppose  that 
Q  6  SO (3 V  I  and  QA  6  svm.  Then  0  =  0o  for  some  line  3  invariant 
under  A  ,  and  in  particular  0  ^  sym. 


P roof .  We  can  suppose  Q  i  I.  There  is  a  unique  :up  to  sirn'. 
unit  vector  x  E  IR°  such  that  Qx  =  x.  Since  QA  S3  sym,  we  have 
QA  *  \QT,  so  QAQ  =  A.  Thus  QAx  =  Ax,  so  Ax  -  cx  for  a 

constant  c.  Each  of  Q  and  \  leaves  V  =  ,  the  orthogonal 

complement  of  x,  invariant,  and  \  is  not  identically  zero  on 

this 

i  i  x 1 


Let 

S 

=  QA  t 

;  s  yn ,  s 

imnl 

ies 

S  ;  V 

=  C.  A  1  V 

i  s 

the 

1  i  ne 

through 

S“  =  A"  .  Since 
g ]  v : n g  0  =  I  or 
x  .  Then  n  ^  sym 


Q’V  is  a  rotation 

^  -  0, ,  ,  whore 

*  i  x  i 

as  in  Lemma  o.l.  ■ 


It  follows  that  if  dim  ker 
QA  £  sym,  then  QA  =  AQ,  so  as 
A  and  Q  can  be  simultaneously 


1  imi 

0  is  both  orthogonal  and  symmetric, 
diagona 1: ted. 


Proof  of  o.l.  If  A  has  J^cinct  eigenvalues,  its  eigenvectors 
are  unique,  so  0  is  either  the  identity  or  is  a  rotation  by  ~  about 
one  of  the  eigenvectors.  M 

Proof  of  0.5.  Suppose  0  a  =  b  *  c ,  and  let  w  be  an  eigenvector 

corresponding  to  the  eigenvalue  c.  Let  V  be  the  plane  orthogonal 

to  >;  V  i  ^  the  eigenspace  with  eigenvalue  a.  As  Q,  A  can  be 

s  :.:)u I *  meoi;  1  -  /  diagonalized  and  Q  is  a  rotation  by  tt  (excluding 

'  =  av  a',  e  either  Q  =  Q  .  or  Q  «  Q.  for  i  a  line  in  V. 

(  W  J  // 

; n  me  farmer  case,  Q.  ^A  has  eigenvalues  (-a,-a,c)  so  has 
:i 1 1  ,i\:s  of  euiii  1  ibriur..  In  the  latter  case ,  Q?A  has  eigenvalues 
a ,  -  i , -c  ;  ,  so  w  is  an  axis  of  equilibrium,  (see  3,3).  B 

o.'.  Corollary.  The  HP*  in  Proposition  6.3  is  a  right  coset  of 

the  subgroup  S*  of  all  rotations  about  w;  in  fact  HP*  = 

w  *  1  0 

where  is  any  fixed  line  in  V,  the  plane  orthogonal  to  w. 

Proof.  HP*  =  (Q  |  V}  and  we  have  the  easily  verified 


identity 


where 


Y) 


Q, 


i/2 


exp*.  ;w 


■*o 


rakes  an  angle 


n  th 


: n  t  he  no  sit: ve  sense  :  n 


These  lemmas  also  enable  us  to  handle  the  next  tvue 


6.8.  Definition .  We  shall  say  A  is  of  type  2  if  A  has  no 
axis  of  equilibrium  and  all  three  of  a,b,c  are  equal  (and  so  ^  o ) . 

6.9.  Proposit ion .  If  A  is  type  2,  then  C\  '  sym  consists  of 
one  point  (A)  and  an  3R1P“, 

Proof.  This  is  immediate  from  6.5.  8 

Notice  that  each  point  of  the  IRP~  has  a  whole  circle  of  axes 
of  equilibria;  namely  Q^A  has  as  axes  of  equilibria  all  vectors 
orthogonal  to  l.  The  eigenvalues  of  Q^A  are  a, -a, -a. 

Types  0,1,  and  2  exhaust  all  symmetric  matrices  with  no  axis 

of  equilibrium,  and  it  is  easy  to  check  from  the  above  that  any 

symmetric  A  with  dim  ker  A  ^  1  lies  on  the  SO (5) -orbit  of  a 
type  0,1,  or  2  matrix.  From  now  on  we  shall  call  these  orbits,  or 
any  representatives  of  them,  type  0,1,  __or  2 . 

Finally  we  turn  to  the  remaining  Afs  with  an  axis  of  equilibr 
that  is  not  already  on  an  orbit  of  type  0,1,  or  2. 

6.10.  Definition.  We  say  A  is  type  5  if  dim  ker  A  =  2  and 

say  A  is  type  A  if  A  =  0. 

6.11.  Propos  it  ion .  I  f_  A  type  > ,  then.  0^  1  >  ^ym  yyny.iy  t  q  __i 
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Proof.  S  =  QA  £  synm  implies  S~  =  A“  and  so  again  S  =  f  \ , 
as  in  0.6  even  though  possibly  A'V  =  0  .  In  this  case  Q  could 
be  any  rotation  about 

All  the  foregoing  information  can  be  summarised  as  follows: 

6.12.  Theorem.  The  SO (5)  orbits  in  M.  fall  into  five  distinct 
types  according  to  the  way  in  which  they  meet  sym  (sec  Table  1  below.)  . 
Furthermore,  for  A  6  sym,  S  ^  =  \QjQA  €  svm}  consists  of  1  {Qi?} 
for  all  l  invariant  under  A  (and  hence  ■“  sym)  except  in  the 
cases  (1)  dim  ker  A  =  2;  then  also  contains  the  rotations 

through  any  angle  about  the  eigen- axis  of  A  corresponding  to  the 
non -zero  eigenvalue; 

(2)  A  =  0;  then  S  =  S0(5)  .  See  Table  2  below. 


Remarks 

1.  Table  1  highlights  the  fact  that  having  an  axis  of  equilibrium 

or  not  is  not  an  invariant  of  the  SO (3)  action  on  L  .  This  means 
there  are  equilibrated  loads  having  an  axis  of  equilibrium,  but 
which, when  rotated  globally  bv  a  certain  amount  to  another 
equilibrated  load A no  longer  have  one. 

2.  Thus,  bi^  Theorem  5.1,  we  get  existence  of  solutions  to  the  traction 
problem  for  all  types  of  astatic  loads  except  3,4.* 

3.  The  notion  of  tvpe  can  he  nulled  back  from  to  L  with  a 

little  care,  as  we  see  below. 


*In  particular  the  occurrence  of  0  solutions  by  Stoppelli  in  *yne  1 
is  seen  to  be  due  to  a  neglect  of  the  full  rotation  group  < see 
Section  S)  .  Our  results  arc  also  consistent  with  those  of  Rail 


6.15.  Definition. 


By  analogy  with  oar  definition 

SA  =  {Q  s  SO (3)  !  QA  e  sym/  ,  A  €  M 
which  we  applied  when  A  is  of  type  0  ,  let  us  now  write 
S(?  =  {Q  e  SO (5)  ;  0^  e  Lo;  ,  l  €  L  . 

From  equi variance  of  k  we  clearly  have: 

6.14.  Lemma .  ^  . 

Note  that  the  map  -►  0^0  sym:  0^  QA  is  an  embedding  for 

types  0,1,2  but  not  all  types  3,4,  because  of  the  isotropy. 

Pulling  back  to  L  ,  we  see  that  Qt—  0£  is  an  embedding  S,  -*■ 

if  k  (Z)  is  of  type  0,1  or  2,  so  we  can  refer  to  <  as  being  of 

type  0,1,  or  2  according  as  k('l)  is.  On  the  other  hand,  if  k  ( ?. ) 

is  of  type  5  then 

e  i  t  her  (a )  0 ,  0  i  =  f  l  y -  l } 

-  /.  e 

or  (b)  0 y  0  L  *  =  two  disjoint  circles  in  *  ker  k 

Finally,  if  k(?,)  is  of  type  4  then  C\  C  ker  k  -  L  and  any 

SO ( 3 1  orbit  C ,  is  all owab 1 e . 


lire  Set  Ligenvalue 


hit  of  A  with  sym. 


Type  of  A  j  Description  of 

i 

four  points 

two  points  and  IRIP*  36  S 

!  one  point  and  !P1P“ 

i 

!  two  disjoint  circles 

!  SO  (3) 

TABLE  : 


Figure  3  illustrates  some  simple  examples  of  loads  of  different 
types.  These  loads  are  all  pure  traction,  with  B  -  0. 


0 

1 

3 

4 


axis  or  equilibrium 


Type  1.  Rotation  by  180° 
about  one  of  the  horizontal 
axes  produces  an  equilib¬ 
rated  load  with  no  axis  of 
equilibrium. 


Type  2 .  Any  horizontal  axis 
is  an  axis  of  equilibrium; 
vertical  axis  is  not  an  axis 
of  equilibrium.  Rotation  by 
180°  about  the  vertical  axis 
gives  an  equilibrated  load 
with  no  axis  of  equilibrium. 


Type  3  (cl)  .  The  load 
itself  admits  a  cireU 
group  of  summe tries 
about  the  axis  t  -  - 
which  is  thus  an  axis 
of  equilibrium. 


Type  3  (  o)  .  The  load  is 
not  symmetric,  but  the 
astatic  load  remains 
constant  under  rotation 
about  the  axis  — which 

is  thus  an  axis  of 
equilibrium . 


I 


i 


Type  .  The  astatic 
load  is  zero:  all  axes 
are  axes  of  equilibrium. 


,oad  types 


3 7 .  The  Bifurcation  Equation  and  its  Gradient  Character 

According  to  the  formulation  (PV,  of  our  problem,  we  wish  to 
solve  the  equation  H(\,A,n;Q)  =  0  for  Q,  where 


Ht\\  , A, n;Q)  =  SkewiQA)  -  \F ( \  , SvnfQAl  , An )  t 


(A,n)  is  near  <A^,n^)  <£  and  A  is  small.  Define  the  right  -  invar:  ant 

vector  field  X.  on  SO (3)  by 
X0 


X  (Q)  =  Skew(QA  ).Q; 


i.e.  right  translation  of  Skew(QA^)  E  so(3)  =  TjS0(31  to  TqS0(3)  . 
Likewise,  we  shall  regard  H  as  a  right- invariant  vector  field  ci 
SO(3)  depending  on  the  parameters  A,A,n  by  setting 


X(\,A,n;Q)  =  H(\  #A,n;Q)  •«} 


Thus , 


XfO,An,n0;Q)  =  x  IQ) 
0 


finally,  note  that  3  is  the  zero  set  of  \  ;  i.< 

A0  A  > 


S  ^  =  *0  SO  (31  Skew  '  v\, , 


What  *3.  is  for  various  type0'  w.a 
\) 


/no  ?  ^  j 
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1 .  Lemma.  Suppose 


a  -  svm 
0 


is  of  anv  tvoe.  Then  por  0-5 


T  S 

q  a. 


{ WQ | W  €  skew  and  WQA0  +  QAqW  =  0;=  ker  DX^  (Q). 


Proof.  The  second  equalitv  is  clear  for  anv  Aft ,  because  DX^  (Q) : 

'  '  U  ;Kn 

*VQ i — *  skew  i.tVQA  )  -Q  .  For  the  first  one,  the  inclusion  C  immediately 

follows  by  differentiation  of  X  (Q)  =  0  in  Q.  Equality  then 

A0 

follows  by  a  dimension  count;  recall  from  3.3  that  v  v  gives 

an  isomorphism  from  the  space  of  axes  of  equilibrium  for  A  (not 

necessarily  of  unit  length)  to  the  W  €  skew  such  that  WA  +  A IV  =  0.  ■ 

T 

Recall  that  W  h-  WQA^  *  W  corresponds  to  the  linear  trans¬ 
formation  trace (QA^) I-QA^  under  the  isomorphism  of  skew  =  S0(3) 

with  IK0.  When  Q  £  S  ,  QA^  is  symmetric,  so  this  transformation 

A0 

i s  sy mme trie  re  1  a t i ve  to  the  Killing  form  on  SO ( 3 ) .  Th i s  rema rk 
and  ".1  proves  the  next  lemma. 
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Lemma.  Suooose 


of  S,  ,  the  range  of 

\) 

complement  of  T . 


is  of  my  type.  Then  at  each  point  0 
DXX  .Q):  T^S0(3i  -  T^S0‘3)  is  the  orthogonal 


Next  we  recall  a  general  context  for  the  bifurcation  of  vector 

fields  that  will  be  applied  to  our  situation  \cf.  Reeken  '  1 9 3  ] j  -  Lot 

M  and  A  be  manifolds  and  X;  M  <  A  -  TM  a  smooth  vector  field  on 

M  depending  on  the  parameters  3  G  A.  We  seek  the  zeros  of  X.  For 

'  =  (y  suppose  the  zero  set  S  of  X  is  a  known  smooth  compact 

submanifold  of  M.  Assume  that  M  carries  a  Riemannian  metric  and 

that  for  x  G  S,  the  range  of  D  X(x,.\q)  is  the  orthogonal  complement 

of  T^S .  The  normal  bundle  E  of  S  trivializes  a  neighborhood  U 

of  S.  For  each  x  G  U,  let  P  :  T  M  -*•  T  S  ,  ,  be  the  orthogonal 

x  x  x  t(x)  b 

projection  to  the  fiber  over  t(x)  ,  where  t:  E  ^  S  is  the 

projection.  By  the  inverse  function  theorem,  there  is  a  unique 
section  :  S  —  E  such  that  P  X(-J>A  (x),\)  =  0  for  x  G  5  and 
in  a  neighborhood  of  (assume,  for  example,  that  ,  *  is  compact  ). 

Let  X(x,A)  he  the  orthogonal  projection  of  X  (x , ;  )  onto  the  tangent 
space  to  the  graph  of  at  a  point  x  on  the  graph.  Thus,  X  f  x ,  '*  1 
is  a  vector  field  on  the  graph  of  f.  and  finding  its  zeros  is  cl  carl}' 
equivalent  ifor  small  '•)  to  finding  zeros  of  X.  We  call  the  equation 
X ( x  i  =  0  on  the  graph  of  t-  the  bifurcation  equation.  Since  ^ 
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and  the  graph  of  are  di f feomorphio  under  ,  we  can  equally 

well  regard  X  as  a  vector  field  on  S.  This  reduction  of  the  problem 

is  often  known  as  the  Liapunov-Schmidt  method. 

The  above  procedure  may  be  applied  to  our  vector  field 

XL\,A,n;Q)  with  parameters  (\,A,n)  and  variable  x  =  Q  <=  SO (51  =  M. 

Assume  \  is  near  zero  and  (A,n)  is  near  a  load  where  A, 

is  of  arbitrary  type.  Thus,  there  is  a  unique  section  ^  of  the 

normal  bundle  to  S.  determined  by  the  Liapunov-Schmidt  procedure 
A0 

as  described  above.  Let  7(.\,A,n)  denote  the  graph  of 

and  let  X(\,A,n;Q)  be  the  orthogonal  projection  of  X  to  the 

tangent  space  of  7  at  Q.  Thus,  X  is  a  vector  field  on  7. 

As  above,  we  may  also  regard  X  as  a  vector  field  on  $.  . 

A0 

The  rest  of  this  section  is  devoted  to  proving  that  the  essential 
part  of  X  is  a  gradient.  In  the  general  context  above,  if  X  is 
a  gradient,  then  so  is  X  since  the  orthogonal  projection  of  a  gradient 
vector  field  to  a  submanifold  is  the  gradient  of  the  restriction.  This 
simple  version  does  not  directly  apply  to  our  situation  as  X  need 
not  be  a  gradient  vector  field  on  S0(5).  However,  the  “second  order" 
Taylor  approximation  of  X  will  be. 

To  state  our  gradient  results,  recall  that  in  ' 4  wo  defined  the 
quadratic  function  G:  L  skew  to  be  the  second  order  term  in  the 
Tavlor  expansion  of  F  about  0.  Thus  F(\,/)  -  <’ !  +  —Cf.M  +  ... 

*  ‘  w  n 

where  G  is  a  quadratic  function  of  l .  The  appropriate  second  order 
approximation  to  the  vector  field  X  will  thus  be  defined  by 

X„r\  , A.n:'}!  =  'skewiOA*  -  -skcw  GiO;n,  •> 


I 


Let  be  the  second  order  approximation  of  the  vector  field  X  on 

S  obtained  bv  the  Liapunov-Schmidt  nrocedure.  Thus ,  X-,fO)  is  the 
*  0 

orthogonal  projection  of  onto  the  tangent  space  TnS  *ov  Q  ^  S*  • 

'  A0  0 


'.5.  Theorem,  Suppose  is  of  arbitrary  type.  Then  X-,  is  a 


gradient  vector  field  on  S  .In  fact,  X_  =  -grad  f,  where 
^  \  _ _ 7  “  7  _ 

*0 


tCQ)  -  Uq.Q  I8>  *<t0,  *W 


<J0'qTi8>  +I 


<7u0’c(7u0)>d' 


and  urt  =  D$(I  )  ;  i.e.,  is  the  uniaue  solution  in  C 

-  Q  v  &  ^0  -  Q  - * -  svm 

of  the  linearized  equations  with  load  Q*,-.  €  L 

- . -  o  e 


Recall  that  the  pairing  between  loads 
(or  displacements)  is  given  by 


;  =  ■  B . 


and  conf igurat ions 


^  B(X)  •  i  (X)dV 

>  8 


'-8 


~  \  ;  • ;  \  dV  =  t  race  k  ( l ,  * ) 


and  physically  represents  a  potential  for  the  working  of  the  loads. 

T  t 

Observe  that  if  l  £  L  ,  then  (  *.  ,0‘l.v  -  trace  '  AO  \  -  tracefAQ  )  =  <  ;  Ql  ' 

e  2  o 

for  al 1  Q  -  SO (5) . 

Remark .  In  the  second  term  of  and  f  we  can  replace  -1 

bv  l  .  However,  the  difference  is  higher  order,  so  is  sufficient 

for  subsequent  applications. 


Proof  of  7 . 3 .  We  shall  show  that  X,  is  a  gradient  field  on 
SO{31  which,  hv  the  remarks  follwoinc  is  sufficient. 

Wc  proceed  in  two  parts.  Let  us  first  show  that  X  '01  is  the 
gradient  of  (  on  all  of  SO (31  . 


7.4.  Lemma.  Let  i  ^  L  and  A  =  kUJ.  Define  i  vector  field 


on  S0(3) 
of  SO (5) 

Proof. 
If  B, 

and 

If  B 


To  prove 


by  X^(Q"i  =  Skew(QA)*Q  as  above  and  the  map 
to  IR  bv  (C,}  =  <  ij  ,  Then  =  -grad 

Two  simple,  but  useful  observations  are: 

,  W  G  M- ,  with  W  G  skew,  then 
o 

(  B,W>  =  (  skew  B , W  > 

G  M  ,  l  G  L  and  t  G  C  »  then 

<  Z,B <t>>  =  <  B , k  ( x  ,  p )  > 

4,  we  compute  as  follows: 

dt(Q)  •  (WQ)  =  <  ?.,  (NQ)TIg  > 

=  <  Ci-.Qj T,  k(.J„)>  bv  (: 

r  o'  ' 

=  <  (WQ)T,A  ) 

=  (  K  r , OA 

=  W ,  skew  i  0 A  ))  bv  »  I 
=  WO,  skew  '0\)*O  • 

=  -'WO,  X  -O;.,  T 


This  deals  with  the  first  term  of  X-,.  To  deal  with  the  second 


term,  we  need  a  special  case  of  Betti's  reciprocity  theorem: 

7.5.  Lemma.  =  <  (WO)Z^fu^>  for  and  fWQ)  ■=  £vm  . 

This  is  a  trivial  consequence  of  svmmetrv  of  DOfl^j  i.e.  of  the 

o 

elasticity  tensor.  It  is  also  proved  in  standard  references;  for 
example,  see  Truesdell  and  Noll  ~ 1965] ;  p.  525. 

To  prove  7,5.,  we  shall  also  need  to  calculate  the  second  derivative  of  the 
skew  component  of  P;  i.e.  of  F(6)  =  Skew [ k ($ (o))j  .  Surprisingly, 
this  second  derivative  depends  only  on  the  classical  elasticity 


tensor 

c  .  Recall  from 

§2  that 

we  regard  c  as 

a  linear  map  of 

svm  to 

itself  and 

that  we 

write  e 

=  L(?U+(7u)T) . 

7.6. 

Lemma.  Let 

F:  C  - 

skew  be 

defined  by'  F($) 

=  Skew  ;k ($(*)) 

Then 

F(TJ  =  0, 

o 

DF(Ig) 

=  0  and 

D“F(I  )(u,ul  =  2  Skew  (!  7u-efe)dV  =  -2  Skew  k(7.  .u) 
b  JB  u 

where  ^  =  -DIV(c(e))  and  7  =  c(e)*N  ,  Identifying 

skew  with  this  becomes 


-D“F(Ig) (u,u) 


b  v  u-  dV  +  f  x  '  n  dA 
u  J  u 
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6 


IS 


Proof.  By  Lemma  4.2,  H  p)  =  Skew f  PdV  where  P  is  the 
'  '  -  B 

first  Piola-Ki rchhof f  stress  tensor.  We  have  Pfl’j  =  i,  so 

S' 

f  -Np  " 

F(I_)  =  0.  Also,  DF(In)*u  =  skew  ,  -—--Tu  dV  -  skew  c*e  d\  - 

S  8  igjF  ;R 

3P  0  2^ 

as  ce  is  svnunetric  and  since  — (  I  )  =  c.  To  compute  o  r,  *e 

j  r  O 

shall  need  to  use  the  fact  that  S  is  symmetric.  Write  P  =  FN  uu: 

use  the  product  rule  to  obtain  DpP(F)»7u  =  7u*8(  Pl  +  FP^Sl'F;  •  .*u. 

Thus,  as  S(T  _)  =  0, 

o 


D“P  !' I_)  •  (7u ,7v)  *  7u-D_SU«)*^v  *  7vDcS(F)*7u  *  DlSHo)  •  ( 7u.7v\ 
8  F  o'  F  F  b 


T  "* 

Now  D-S(I-)#7u  =  DrS(I„)  •  (7u  +  7u  )  =  c*e  and  D  S  {  I  n)  is  symmetric, 
r  o  C  o  r  o 

SO 


D:FfIJ-(u,u) 

o 


skew 


!  DjIT  ( I R)  ( 7u ,  7U)  dV 

j  r  o 

8 


=  2  skew 


(I  7u#c(e)  dV 

•j8 


Finally,  this  equals 


-  2  skew  { j  bu®  u.  dV 

j 

B 


j  vc  dA  } 

^8 


by  the  divergence  theorem,  so  the  last  statement  follows,  f 


50 


-GU)  =  2  Skew  [  j  b  ®  u^dV  +  j  r  ®  u?dAj 
S  cB 


=  2  Skew  k((b,x) ,u^) 


where  b  =  -DIV(c*  (e)  J ,  and  t  =  c(e>)*N  and  =  [7u?  r  (7u;,  )  }  ■' 2  . 
However,  these  last  equations  say  exactly  that  (b.i)  =  l ,  and  so  we  get 


jGU)  =  Skew  kU,u  ) 


Completion  of  the  proof  of  7.3,  The  derivative  of 
\  T 

Q  >*•  (  VQ  UQ>  in  t^ie  direction  WQ  is  given  by 


<  st0.,  y(wQ)TuQ>  *  <  i.Q,  W^m  > 


<  ,  (WO)  u^)  (by  Betti  reciprocity,  ~.51) 


( V"o 


.<  V  ,  k  (Q?-0  .  uQ)  > 


bv  (2) 


■  <  V!  ,  Skew  k  ( ,  u^)  )  by  (1) 


-<  WO  ,  Skew  k(QcQ,Ug)Q  ) 


<  tvq,  ±OfQ£  5-0  > 


by  i  3 )  M 


38 .  Bifurcation  Analysis  for  Type  I 

We  now  discuss  the  solutions  of  the  basic  equation 

H(a , AfM;0)  =  Skew (OA)  -  \F(a,S>tti(OA]  ,Qn)  =  0  i. 1  ) 

for  the  load  i  =  (A,n)  near  a  load  having  an  axis  of 

equilibrium  and  of  type  1,  and  for  A  near  0  .  We  shall  also  obtain 

the  stability  of  the  solutions  and  finally  we  shall  compare  our  results 

with  those  of  Stonpelli  [1958].  For  type  1  we  need  to  do  a  bifurcation 

analysis  on  the  circle  S.  corresponding  to  the  degenerate  zero  set 

A0 

of  H  when  a  =  0  and  l  =  .  The  analysis  has  some  features  in 

common  with  the  papers  of  Hale  [1977]  and  Hale  and  Taboas  [ J  9S 1 ] . 

Without  loss  of  generality  we  can  assume  that  A^  =  dia? r a , -  a , -c ) 

'  'y 

where  0  f  a~  f  z~  *  Thus,  from  §6  the  set  3.  of  zeros  of  Skew  (OA)  for 

A0 

0  6  S0(3)  is  given  explicitly  by  the  following  two  points  and  circle: 


S  -  [ diagf 1  *  - 1 , - 1)  diag(-l , 1 1) }  U  C  (2) 

An  "  Ao 

where 

x  -y  0  \ 

y  x  0  J  'x  =  cos  A,  v  =■  sin 
0  0  1  / 


The  loads  corresponding  to  the  two  points  are  A  =  diag  <a,a,c‘.  and 

+■  * 

A  =  d iag i -a , -a , c )  . 

From  ".5,  we  are  led  to  study  the  critical  points  of 

T  T 

f<Q)  =  (  <  ,0  )  on  .  Note  trom  the  vlivergencc 


theorem 


hat 


Thus  the 


<  lo*\  >  = 


•un)  c(ej  >dV 


where  =  D$(I)  (0?-q)  and  =  [  >u^  *  ('Tu^j‘1/2  .  Thus  the 

function  f  is  computable  from  linearized  elasticity  alone,  which 
leads  to  the  curious  observation  that  our  ’’second  order'1  nonlinear 
elasticity  here  involves  no  more  data  than  linear  elasticity,  but 


merely  processes  the  information  in  a  different  way.  Writing  Q  = 
/  x  -y  0  \ 

f  y  x  0  J  as  in  (2)  ,  f  becomes  a  nolvnomial  of  degree  2  in 

\  0  0  1  J 


(x,y).  Write  the  two  terms  of  f  as 


f  (Q)  =  tU,y)  = 


+  blx  ♦  b:y)  +  ^arx- 


a7xy  *  ay-  ♦  ax  *  a.y  +  a  ) 


which  defines  the  numerical  constants  by.bj.b.,  and  al . ao‘  Sext 

define  new  parameters  ot^  by  writing 


f  (x,y)  =  ff(x,y) 


and  letting 


f  (x,v)  =  ijx"  ♦  a-jxy  ♦  u..y"  +  XjX  ♦  a.y  +  :i( 


Note  that  x  . . . ,  x^  depend  on  our  parameters 
the  elastic  moduli  of  the  material.  Thus, 


as  well  as  on 


‘i  =  V  a:  =  a:’  ‘5  =  a:> 


=  ^b,  +  a,»  i-  =  +  :u,  ’  V,  ~  *  ‘h. 


L1  "  •  1  V  3 


Replacing  •  ^  by  ,  where  Q  is  :xs  in  '2),  effects  a 

rotation  of  the  x-y  plane.  Thus,  by  rotation  of  r  if  necessary, 
we  can  assume  =  0  . 

Let  us  fix  i.  ,  x,  and  consider  the  bifurcations  of  zeros  of 

1  j 


df* 

db 

-  2( } xv  + 
o  1 

x_x 

-  1  y  on  ('i.e.  , 

cri t leal 

points 

of  f* 

on 

i 

5  * )  with  x . 

4 

and 

i-  as  parameters. 

Set  M  =  1 , 

,  zu  ,  b 

D 

>  c  ci  <  df*  r 

'  dr  4 

=  o; 

•  ,  the 

mani fo  Id 

of  c 

ritical  points 

of 

f * .  Indeed,  M  is  a 

man i f o 1  \ 

and  can 

be 

parametrized  by  o 

:  IR 

<  s1  -  m  ,  = 

(-Zfzij+u) 

cos  f , 

~  ^  "3 

2  12 

Denote  by  t:  IR  *  S  •+  1R  the  projection  onto  the  first  tactor. 


3.1.  Lemma.  Set  2  -  f2(a  -ou)"  -  a  -  108a  i  -:0“  .  ("*) 

-  -  1  o  a  d  4  ;>  1  :>' 

If  a,  -  'a.  ^  0  ,  then  t.  M  -*■  is  a  proner  stable  nan  in  f t ,, 1  -  space  , 
—  1  j  -  - 4 - * - * -  4  o  — - - - 

and  its  set  of  critical  values  is  the  astro  id  defined  by  =  () 
i  see  Figure  **.  bel ow  )  . 


Since  the  number  of  points  in  ^  ft}  (i.e.,  the  zeros  of 


df* 
d  ■ 


at 


is  a  constant  over 

4  o 


0 


or 


0  , 


we 


obtain 


S.2.  Corollary. 


df* 


has  4  zeros  i  f 


0  ,  and.  has  2  zero-; 


0 


■ 
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Proof  of  Lemma  3.1. 


The  critical  set 
A 


of  t^o:  R  <  S 


1 


R" 


i  (a,6)  £  R  <  S  !  ot  sin”9  +  cu  cos“8  +  y  =  0}  . 


critical  values  of  tt  can  be  parametrized  by 


is 

Thus,  the  set  of 


jl  =  -2  (a,  -ou)  cos°  3 
4  1  o 

^5  *  2(ctj-ot,)  sin^  9 


rs) 


Since  I  consists  of  4  cusp  points,  4  fold  lines  and 
voc jl  is  1-1  ,  by  a  result  of  Whitney  (see  Mather  [1969]  or 
'^olubitsky  and  :uillemin  [1973]),  one  knows  that  7Top  is  a  stable  map. 

Eliminating  9  produces  the  bifurcation  set 
r^r  N12/3  2/3  2/3 

[2(Dt ,-cu)]  =  a,  *  (9) 

1  O  4  D 


For  -  ou  /  0  ,  (S)  describes  the  astroid  shown  in  Figure  *4. 


Figure  4- 


Next,  observe  that  tor  real  numbers  A,  R  and  C, 


A  +  B  +  C  =  0  if  and  only  if  A*3  +  Rv>  +  (T3  =  3ABC 


by  virtue  of  the  identity  +  R'  ♦  C''*  -  3ARC  =  •  A  ♦  R  «■  ,'ifA" 

AB  -  BC  -  CA)  .  Applying  (101  to  shows  that  :  0  1  i  oquivalon 

2  :  :  :/s  j/3 

a ,  +  'i„  -  2  (a,  -  ju )  =  -3ft,  x.  (-ft,  -  i  - i 

4  ^  l  a  4  .)  1  a 

gives  the  stated  conclusion.  M 


fubmr  both 


The  family 


df* 

d9 


of  functions  on  S 


with  parameters 


A- 

o 


so 


^4  ' 


eniovs  a  universal  property.  Consider  a  perturbed  family 

gf2-  +  gfv.c.e)  ,  With  g(0,0,^  =  0  for  O.c)  ^  IR  ‘  T™  .  To 

,  df* 

each  (\,c),  denote  by  (-p—  +  =  0; 

the  "manifold"  of  zeros. 


S . 5 .  Lemma .  For  (\,c)  sufficiently  small,  the  sets  M.  %  are 

manifolds  and  there  exist  two  smooth  families  of  diffeomorphisms 

:  1R“  1R“  and  ‘f.  :  M,  -  M  defined  for  \,c  sufficiently 

•  ,  C  -  A  ,c  *  ,  C  -  - 

small,  such  that  nV.  =  *>.  ,  and  n  =  identity,  v'  .  =  identity. 

— ■ — — 2 -  a  ,  c  •  X  ,  c  — —  0,0  '0,0 


1  1 

Proof:  For  \,c  sufficiently  small,  the  map  :  IR  *  S  -*  M, 

-  '  'A 

cu  c(y;a)  =  (-2(o^+u)  cos  +  sin  9  ,  -2(ou+u)  sin  ^  -  cos  ) 

defines  a  parametrizat ion  of  M.,  c  .  By  Lemma  S.l,  11  '  S*  -*  F 

is  an  unfolding  of  the  prone r  stable  map  .  Tims,  one  can  find 

1  2 

di f feomorph isms  if*.  ^  ^  on  ]R  <  S‘  and  F“  respect  i  vely 

such  that  *0) 9  f  *,  .  This  lemma  follows  bv  letting 

A,c  A  ,  C  '  ,  C 

?.  =  .v-f*,  ■ ; 1  .  ■ 

’  \  ,  C  *  *  ,  C  A  ,  C 


\'ow,  we  ire  ready  to  state  our  main  result  on  the  number  of 
solutions  near  t  ype  1  loads.  Let  -  =  •'  t  c  •  depend  smoothly  uj  a 

parameter  e  in  F™  ,  with  •  ■. n  1  ~  .  Pe/ail  that  A  :  "• 

by  ■  1  i  . 


le fined  bv 


a,,  a ^  and  a. 


3.4.  Theorem .  Let  se  a  type  1  load  with  k i 4  ;  =  , a , -a ,  -c  )  , 

*»  i 

0  i  a”  c“  ,  a-,  =  0  ,  and  a.  i  .  Then,  there  exists  a 

2  -  1  j  - 

(smooth]  function  A (a, cl,  A(a,0)  =  d.(a^,a„]  +  0(‘*J  defined  for 
(.  \ ,  c  1  sufficiently  small  and  \  >  0  such  that,  the  traction 
problem  has: 

( i )  four  solutions  for  the  load  \  1(c)  ijf  A (  •  , c  !  * '  0 

i two  of  them  near  C  )  . 

A0 

( i  i )  six  solutions  for  the  load  \  Z  (c )  i_f  A  0 ,  c )  >  0 

(two  of  them  near  C  ]  . 

A0 

Proof:  The  bifurcation  of  zeros  of  a  (cf .  3~)  on  C.  (=$*) 

A0 

is  the  same  as  finding  zeros  of  (  -X,  4—  >  =■  —  '  +  —  g(\,cf?) 

d  f  df  * 

=  —  +  ~  g  or  -rr  +  2  ,  where  g(0,0,f]  =  0  .  Let 

dtr.  2  df  &  &  * ,  c 

be  the  family  of  dif feomorphi sms  found  in  Lemma  S.3.  Take 
bj(&(c))  b^(s(c)) 

AO  ,  c  ]  =  A Ou  c  (  — - -  +  a4  ,  ^ -  +  a. ]  which 

has  the  desired  property.  1 


Next,  we  want  to  determine  the  "generic”  structure  of  the  bifurcation 

set  K  =  ;A  =  0;  in  (>.  ,c)  space,  '  >  0. 

F;or  m  =  0  ,  AOc  ,/:u)  =  A(a,,a~)  ±  0  ,  it  is  clear  that  K  -  $  . 

4  o  4o 

Indeed,  our  traction  problem  has  two  solutions  near  C  it 

0 

A  (a  ,a.l  •  n  ,  and  four  solutions  near  C  if  fa1ta^)  0  . 

1  j  \j  • 

b  A  fc  :  :’.Vcl 

For  m  =  1  ,  consider  k:  c  »— *  1 — -  *  a  ,  — ■ 


where  \(c!  is  the  linear  part  of  ; .«:  c  '  .  represents  a  straight 

line  which  we  assume  to  intersect  the  astroid  t  ransv^rse  1  v  if  they 
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meet.  Notice  that  K.  -  •  c  '  <  * . ,  c  ^  <  is  the  inverse  image  of 

\ 

the  astro  id  (defined  bv  equation  (9)),  under  the  man 
brUc)  '  ‘  b  ,l(c) 

hy  :  c»—y,  ( — r-r -  +  a.  ,  — ^ -  *►  a.  1  .  f!j,  as  in  the 

a  ‘,c  24  4  2  a  o  ■  a,c 

proof  of  Theorem  3.4).  Write  \c  =  c  .  Recall  that  ?.  (cj  =  Afcj  +  Of 

and  consider  the  map 

b  A(c)  by\(c) 

h,  :  c  H-h,  (a,c)  =  ,  -  ( — ^ -  +  a  +  Q(X)  >  -=-? -  +  a.  +  0(A))  . 

^  '  ,  ac  —  a  —  d 

Since  the  astroid  is  bounded  and  x  is  close  to  the  identity, 

a  ,  c 

there  exists  an  interval  (-M,M)  such  that  K ,  =  ;c\c  -  <K.  }  C~  f-M,M) 

'  \ 

for  V  >  0  ,  c  sufficiently  small.  Applying  the  isotopv  theorem 
for  transversal  maps  (see  e.g.,  Hirsch  [1976))  to  the  family  h. 
through  hg  =  k  ,  we  conclude  that  the  bifurcation  set  X  consists 
of  0,  2,  or  4  curves  with  slopes  k  ^(astroid)  (see  Figures  £,&)  . 
Thus,  for  example,  by  choosing  c  i  0  sufficiently  small,  and 
letting  \  0  (i.e.,  consider  the  load  a£(c)),  one  can  pass  from 

a  parameter  region  where  there  are  two  solutions  near  the  circle 
(four  in  all)  to  one  where  there  are  four  near  the  circle  (six  in 
all)  and  back  again  to  the  two-solution  region  (see  Figure  t )  . 

Such  a  situation  is  not  dealt  with  in  the  analysis  of  Stoppelli  [1958]. 


For  m 


let  us  suppose  that  the  affine  map: 


b.A(c)  bnA(c) 

i — -  +  a  .  — : ^ -  +  a_)  is  onto,  where  again  Afc) 

2  4  2  d 


is  linear  part  of  ?,(.c)  .  Without  loss  of  generality,  we  may  also 
assume  that  b^A(c)  =  c ^  and  b^A(c)  =0-,  t  where  c  . 

Notice  that  K,  ^  =  {  (c. , c 1  (\  , c, , ,  z)  -  ^  ;  is  the  inverse  image 

'•  t  -  l  «.  i  ~ 

b1 Af *  ) 

of  the  astro  id  under  the  map  h,  :  (c .  ,e*  1  *-*’  {  -*  - - 

\,Z  1  ~  •  ,c^  ,  -.V 


bv\(:0 

Ar —  *  a, 

-  A  D 


Set 


\c,  r  c.  ,  a c n  =  ,  and  consider  the  man 


is  bounded  uniforml 


As  before,  -C  ,  =  *  (c.  ,cj  p-c  ,  \cp  ^  -A . 

\ —  i  -  ^  - 

for  \  >  0  ,  c  sufficiently  small.  Applying  the  isotopy  theorem 

for  transversal  maps  to  the  family  h,  ^  through  the  affine 

isomorohism  hA  -  we  conclude  that  the  bifurcation  set  is  a 

0,0 

cylinder-1  ike  set  along  the  c-axis  with  base  a  cone  over  the 
astroid  in  c^,c.,  space.  The  first  order  approximation  of  this 
cone  is  given  by  the  cone  over  the  astroid  in  the  plane  \  =  1 , 
centered  at  (-2a.,  -2a-)  with  "sice’'  4>a1-a-  (see  Figure  7). 


(1 


Next  we  discuss  the  stability  of  the  solutions  corresponding  to 
loads  near  load  of  type  1.  This  can  be  determined  by  combining 
our  stability  results  for  type  0  : Theorem  5.5;  together  with 
well-known  stability  results  for  the  cusp.  We  make  the  same 
assumptions  as  those  in  Theorem  S.4. 


3.5.  Theorem.  Let  A  =  diag (a , -a , -cl ,  A^  =  diag(a,a,cj  and 
*  ■* 

A0  =  diag(-a, -a,c)  as  above .  The  indices  of  the  bifurcating 

solutions  are  boxed  in  Table  3.  f Recall  that  stable  solutions  have 

index  =0.)  In  each  case  the  circle  represents  C.  defined  by 

0 


equation  (2). 


Note  that  two  stable  solutions  bifurcate  off  the  circle  when 
c  >  'a!  .  In  all  other  cases  the  solutions  near  the  circle  are 


unstable . 


3.4.  Example .  Let  8  C  ]R  be  a  region  with  unit  volume  and  let 
the  load  be  given  by  l  =  (0 , where 

*>  n 

0  ^  a"  ?  c’. 


1 


where 


- ( t  race 


Hence 
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which  is  a  constant  ('independent  of  x,  y) .  In  this  situation, 
=  'x  =  0  and  so  our  theorems  do  not  apply. 


3.5.  Example .  Consider  the  same  traction  problem  as  above,  but 
with  a  homogeneous  nonisotropic  hyperelast ic  material  whose 
linearized  elasticity  tensor  is  given  by  cfc)  =  G  -  t  iiag  *2  . 


7  *1  7  7 

4  a  “  x  ~  +  2awy”  +  2c“ 


1 "  ,  , 
Hence,  A  =  3a  “  0  ,  and  our  traction  problem  in  nas 

six  solutions  ("four  near  C.  )  ,  with  stability  determined  by 

An 


table  3. 


Next  we  shall  discuss  how  to  obtain  the  results  of  stopnelli 
1958  .  as  a  special  case  of  our  analysis,  ive  refer  the  reader  to 
the  statements  of  Stoppelli's  results  in  Grioli  ^19b2,  p.  58,.  In 
this  approach  one  focuses  attention  on  bifurcations  that  occur  on  the 
circle  by  examining  what  happens  at  a  particular  location  on  the 
circle  and  ■  -  «-n  .  We  can  assume  this  point  is  i  1 ,'))  i.e.  :  =  ■) 
w:rh  no  loss  of  genera  1 i t y . 

First  of  all,  if  +  ju  t  0  then  (1,0)  is  not  a  critical 
point  of  f * ,  so  there  are  no  solutions  near  (1,0).  We  may  assume 
then  that  xn  +  ~x„  =  0  and  then  the  Taylor  expansion  of  f*  about 
0=0  becomes 

f*'J>  *  hi  *  V  *  ‘-n  *  H  -  r»^  -  T3 


(higher  order  t  o  r ms 


For  critical  points,  we  are  seeking  zeros  of 


l4  jf*  1  \ 

Case  1.  If  -i.  +  y  .  *  —  i  o,  then  -- —  =  2\-<.  +  + 

-  1  o  J  a»-  1 

and  so  there  is  just  one  solution.  This  is  Theorem  I  on  p. 
of  Grioli  1905  . 


Case  2. 


If  -:u  +  -  —  =  ()  and  a.,  *  IJ, 

1  j> 


jf*  ,  _ 


01' 


and  so  there  are  »  ,  ’ 


;  fold  point).  This  is  rheorem  1  on  p .  3S  :>f  lr:oi; 

^4  li  . 

Case  5.  It  -  '4.  ^  =  d ,  a-,  =  0  but  t,  -  4-  *  —  *  1 

-  1  a  2  2  i  o  s 


then 


dr* 


4  .3 


1 


so  there  are  ; 


)r 


d*3  =  5Utl  '  '43  S  ,r 

solutions  (cusp  point).  This  is  Theorem  fl  on  n.  3S  of  Grioli 


Furthermore,  if  we  express  our  constants  i=a.  1  :n  tor: 
elasticity  tensor  c  and  solutions  of  the  linearized  probiem 
using  (3)  above,  we  find  the  same  conditions  for  those  three 
as  is  given  on  p.  S~  of  Grioli  1 4  h  J  , . 

Thus  we  recover  the  results  of  Stoppelli  on  type  1  loads 
was  explained  in  the  introduction,  however,  this  analysis  is 
local  on  the  circle  and  does  not  give  the  full  story  of  the 
bifurcation  picture,  even  in  this  case.  i'he  complr*. e  bifurca 
analysis,  including  stability. is  summarized  by  our  Figure  ,  i 
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